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Abstract 

We consider single-channel transmission through a double quantum dot system consisting of two 
single dots that are connected by a wire and coupled each to one lead. The system is described 
in the framework of the S matrix theory by using the effective Hamiltonian of the open quantum 
system. It consists of the Hamiltonian of the closed system (without attached leads) and a term 
that accounts for the coupling of the states via the continuum of propagating modes in the leads. 
This model allows to study the physical meaning of branch points in the complex plane. They 
are points of coalesced eigenvalues and separate the two scenarios with avoided level crossings and 
without any crossings in the complex plane. They influence strongly the features of transmission 
through double quantum dots. 
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I. INTRODUCTION 



The phenomenon of avoided level crossing (Landau-Zener effect) is studied theoretically 
as well as experimentally for many years. It is a general property of the discrete states of a 
quantum system the energies of which will never cross when the interaction between them 
is nonvanishing. Their wave functions are exchanged at the critical value of a certain tuning 
parameter where the avoided crossing takes place. The reason for the avoided crossing of 
two discrete levels follows from the expression for the two eigenvalues e± of the Hamiltonian 
of the system, 

ei + 62 , 1 



e± = ^Y^±2\/(ei-e2)^ + 4a;2 

where Ci and 62 are the energies of the non-interacting states and uj is their interaction. 
Since the square root contains only positive values, e+ and e_ are always different from one 
another with the only exception of vanishing interaction uj. 

A crossing point of the two eigenvalues e± can be found by continuing into the complex 
plane, i.e. by adding a negative term into the square root. The mathematical properties of 
such a crossing point in the complex plane are discussed in many papers. According to Kato 
they are called exceptional points, since the spectrum is supposed to be incomplete at 
these points. The exceptional points are branch points in the complex plane . Although 
the number of these points in the complex plane is of measure zero, their meaning for physical 
processes is large. They are related to the phenomenon of avoided crossing of discrete states 
as shown already in 

In recent studies, it turned out that not only discrete states avoid crossing, but also 
resonance states do not cross, as a rule An avoided level crossing in the 

complex plane is accompanied by a redistribution of the spectroscopic properties of the 
resonance states. Most interesting is the bifuraction of decay widths related to the avoided 
crossing of levels in the complex plane since it causes different time scales in the system. 
The long-lived (trapped) resonance states cause narrow resonances in the cross section on 
a weakly energy dependent background induced by the short-lived resonance states. A 
similar situation is discussed recently in . The resonance trapping phenomenon discussed 
is a collective coherent resonance phenomenon as stated also in The 
avoided level crossings may form a branch cut This cut can be traced up to the avoided 
crossing of discrete states ^j. 



Often, the branch points in the complex plane are identified with double poles of the 
5* matrix ^ when related to physical processes. It became possible directly to study the 
spectra of atoms in the very neighborhood of double poles of the S matrix by means of laser 
fields 

mm 

. The results show a smooth behavior of the observables when crossine; 
the double pole by tuning the parameters of the laser field. Moreover, recent studies in the 




ane 



framework of schematical models have shown that the branch points in the complex 
separate the scenario with avoided level crossing from that without any crossing 
In 

mm 

, the double poles of the S matrix are identified with points at which the 
eigenvalues of two states of the effective Hamiltonian coalesce. In 

mm 

, the line 

shape of resonances in the very neighborhood of double poles of the S matrix is studied. 

In the S matrix theory is applied to the transmission through double quantum dots 
(QDs) consisting of two single QDs and a wire connecting them. The study of these artificial 
molecules is of great interest since they display the simplest structures of quantum-computing 
devices that can be controlled by external parameters, e.g. jl?!. Iisj]. One of the interesting 
results obtained for a double QD system, is the appearance of transmission zeros of different 
order at energies that are related to the eigenvalues of the Hamiltonians of the single QDs 
[igI i. They appear even in cases when the transmission is large in this energy region. In such 
a case, they can be seen as narrow dips in the transmission probability. 

Double dot systems provide a very powerful tool for studying the properties of branch 
points in the complex plane and their physical meaning. When leads are attached to them, 
the double dot systems allow further to study the relation of the branch points in the complex 
plane to both the double poles of the S matrix and the points where two eigenvalues of the 
Hamiltonian of the open quantum system coalesce. This is, above all, due to the symmetries 
involved in the system in a natural manner. Moreover, the properties of a double dot system 
can be controlled by external parameters in a very clear manner. The double QD itself is 
characterized by the coupling strengths u between the wire and the single QDs, the spectral 
properties of the two single QDs, as well as by the length and the width of the wire. The 
coupling of the double dot system to the environment is given by the coupling strength v to 
the leads attached to it. All these parameters are well defined and can be controlled. One 
may call v the external coupling of the double QD system (via the leads) and u the internal 
coupling (via the wire) that is characteristic of the double dot system as a whole. 

In the present paper, we will study a simple model for a double QD system with the 
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aim to receive information on the branch points in the complex plane and their relation to 
physical processes such as transmission. We use S matrix theory combined with the method 
of the effective Hamiltonian which consists of two parts. The first part is the (Hermitian) 
Hamiltonian of the closed system and the second part is an additional (non-Hermitian) 
term that takes into account the coupling of the states of the system via the continuum. 
The continuum is given by the modes propagating in the two half-infinite Id-leads when 
attached to the system. The interplay between these two parts of the effective Hamiltonian 
characterizes the different physical situations. 

In Sect. II, we give the S matrix for the transmission through a model double QD system 
by using the effective Hamiltonian formalism. The double QD consists of two single QDs with 
one state in each, a wire with a single eigenenergy that depends on the length of the wire, 
and with one channel for the propagation of the mode in the attached leads. We define the 
spectroscopic values and of the resonance states k. In Section III, we study analytically 
the features of the eigenvalues and eigenvectors at the branch point in the complex plane. 
Here, at a certain energy E — Ec, two eigenvalues of the effective Hamiltonian coalesce. 
We show numerical examples obtained for branch points in the complex plane as well as for 
the transmission through the double dot system. The branch points can be seen by varying 
different parameters. The transmission scenario at small v/uis characterized by transmission 
peaks which are spread over a certain energy region that is the larger the larger the internal 
interaction u is. In contrast to this picture, the transmission peaks are no longer spread 
in energy when v/u is large. Here, level attraction and width bifuraction take place with 
the consequence that one narrow resonance appears on a smooth background created by the 
two broad resonance states. The separation between the two different scenarios is provided 
by the branch point in the complex plane. This separation is independently of whether the 
eigenstates cross or avoid crossing in the complex plane at the energy E^. 

In Sect. IV, we consider the effective Hamiltonian as well as the transmission through 
the double dot system when it is coupled with different strength to the two leads. In the 
following section V, we show numerical examples for transmission and eigenvalue trajectories 
of a double dot system with altogether five and eleven, respectively, states as a function of 
both, the length of the wire and the (external) coupling strength v. The main features of 
the eigenvalue trajectories as well as of the transmission are the same as those discussed 
in Sect. HI. Moreover, we draw some conclusions on the different bonds of the two single 
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QDs in the artificial molecule. The appearance of different bond types is also related to the 
positions of the branch points in the complex plane. In the last section, we summarize the 
results obtained. 



II. EFFECTIVE HAMILTONIAN AND S MARIX THEORY FOR TRANSMIS- 
SION THROUGH COUPLED QUANTUM DOTS 

In our study, we follow the paper where the S matrix theory for transmission through 
QDs in the tight-binding approach is formulated, and the paper ^ where the S matrix 
theory is applied to a double QD system consisting of two single QDs coupled to each other 
by a wire. As in we consider a simple model with a small number of states in each 
single QD and one mode propagating through the wire. This simple model is able to explain 
the characteristic features of the transmission through realistic double dot systems of the 
same structure, as shown in jioj l. 

First we will consider the simplest case with only one state ei in each single dot and one 
mode e(L) propagating in the wire of length L. The wire and the single QDs are coupled 



C=L,R 



by u. The effective Hamiltonian of such a system is 

HeS = Hb + ^ Vbc 
c 

where 

He 



16 



3 



E+-Hc 



-V( 



(1) 



V 



El U 

u e(L) u (2) 

U El ^ 

is the Hamiltonian of the closed double dot system. He is the Hamiltonian of the left (C = L) 
and right (C = R) reservoir and = E+iO. The second term of ifgfr takes into account the 
coupling of the eigenstates of Hb via the reservoirs when the system is opened. It introduces 
correlations between the states of an open quantum system that appear additionally to those 
of the closed system [3]. The effective Hamiltonian H^s is non-Hermitian. 

The coupling matrix between the closed double dot system and the reservoirs can be 
found if both are specified. We take the reservoirs (leads) as semi infinite one-dimensional 

n 

wires in tight- binding approach [19||. The connection points of the coupling between the 
system and the reservoirs are at the edges of the one-dimensional leads. Then the coupling 
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matrix elements take the following form |16l . [1 



sin k 



Vm{E, L) = vi)E,L{xL)i)m{j = I) = V]l 



sin k 



Vm{E,R) = VlljE,L{^R)'ipm{j = 3) = V\ — 1pm{3) 

V 2n 



(3) 



where k is the wave vector related to the energy hj E = —2 cos A;, 4'm{j),j = 1,2,3, are 
the eigenf unctions of and v is the hopping matrix element between the edge of the lead 
and the QD. The v will be varied in our calculations. The eigenvalues of the Hamiltonian 
(121) are real. 



Ei,3 - ^ TV, ^2 

with f]^ = Ae"^ + 2u^, Ae 



ei - t{L) 



(4) 
(5) 



and the eigenstates read 

/ 

1 

|1) = 



2?7(?7 + Ae) 



—u 
rj + Ae 
—u 



|2) 



/ 1 \ 




/ u \ 



|3) 



1 



2r/(77 - Ae) 



T] — Ae 



V 



u 



(6) 



As a result, we get the following expression for the effective Hamiltonian 



1 n(n-^. 



H, 



eS 



r?(»?+Ae) 



ei- v'^e'-^ 



V2v 




v 



V2rj 



-^3 viv-M I 



which is symmetric. Its complex eigenvalues and eigenvectors \k) are 



(7) 



Z2 = ex- v^e^, 

ei + e(L) - v'^(^^ 
zx-, = T 



\ 



e{L) — ei + v'^e''^ 



2^2 



(8) 



and 







1 



|3) 



/ h \ 


\-aj 



(9) 
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where 

a = , , = W — 

/=^. . = -„+^^. ^^ = .= + f. (10) 

The eigenf unctions are biorthogonal, i^effl^) = Zk\k) with [3| 

(fc|/) = (r|/) = 4,;. (11) 

Using the eigenvalues (jH} and eigenfunctions Q of the effective Hamiltonian, the amphtude 
for the transmission through the double QD takes the simple form 

Substituting Q, (0) and (jHl) into the matrix elements (L|V^|A) and (A|\^|i?) we obtain 



{L\V\2) = ;^(E,L|V|m)(m|2) = ^J^, 



i2\V\R) = Y.mm){m\V\E,L) = ~lf-^, 



{L\V\1) = il\V\R) = vf-^iMl)a + Mm, 



{L\V\3) = {3\V\R) = v^'-^{Ml)b-Ul)a). (13) 

The transmission probability is T = |tp. 

The spectroscopic values such as the positions in energy of states are originally defined 
for the discrete eigenstates of Hermitian Hamilton operators that describe closed quantum 
systems. The decay widths do not appear explicitely in this formalism since the eigenvalues 
of the Hamiltonian are real. They are calculated from the tunneling matrix elements by 
means of the eigenfunctions of this Hamiltonian. The corresponding values for resonance 
states are energy dependent functions since the eigenvalues as well as the eigenfunctions of 
the non-Hermitian effective Hamilton operator (0) depend on energy. Nevertheless, spec- 
troscopic values for resonance states can be defined, and that by solving the fixed-point 
equations ^] 

Ek = Re{zk)\j,^j,^ (14) 
7 



and defining 



Tk = 2Im{zk)\E^E^ . (15) 

The values Ek and Tk cliaracterize a resonance state wliose position in energy is Ek and 
whose decay width is Tk- This resonance state causes a resonance of Breit-Wigner type in 
the cross section when it is well separated from other resonance states. In the regime of 
overlapping resonances, the relation between Ek and Tk on the one hand, and the resonances 
seen in the cross section on the other hand, is less well defined. 

In the denominator of the S matrix, the eigenvalues Zk of the effective Hamiltonian i^eff 
appear in their full energy dependence. That means that at every energy E of the system, 
the contribution of every resonance state k is taken into account in correspondence to the 
value Zk{E). This fact becomes important when Zk{Ei^k) 7^ Zk{Ek) and the contribution of 
the resonance state k can not be neglected at the energy E = Ei, i.e. when the resonance 
states overlap. 

Another definition of the spectroscopic values of a resonance state is by means of the poles 
of the S matrix. This (standard) definition of the spectroscopic values in the framework of 
the 5* matrix theory is not a direct one since the poles of the 5* matrix give information on the 
resonances, but not on the spectroscopic properties of the resonance states. The S matrix 
has a pole only when the energy is continued into the complex plane. We remind however 
that the S matrix describing physical processes is defined for real energies E, and 15*^ < 1. 
It is not surprisingly therefore that the two definitions do not coincide completely. In the 
following, we will characterize the resonance states by the energy dependent eigenvalues Zk 
and eigenfunctions \k) of the effective Hamiltonian H^s as well as by the values Ek and Tk, 
but not by the poles of the 5* matrix. The reason for doing this, is the clear definition of the 



3, 



by means 



spectroscopic values Ek and Tk also in the regime of overlapping resonances 
of the effective Hamiltonian i^efr that describes the open quantum system. 

It may happen that, at a certain point, Zk = zi for two different states k and I. Such a 
point might be considered as the analogue of a double pole of the S matrix. However, the 
coalescence of two eigenvalues Zk, zi at a certain energy Ec does not mean that also the 
poles exactly coincide. Therefore, we will not consider double poles of the S matrix in the 
following, but will look at the points and their energies Ec where the two eigenvalues Zk, zi 
coalesce. In such a case, the transmission is determined mainly by interferences between 
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the two resonance states k and /. These interferences influence strongly the line shape of 
resonances 0, [3| ■ 

Generally, two resonance states k and / avoid crossing in the complex plane, i.e. the 
eigenvalues Zk and zi coalesce at an energy E = that is different from the energies i^^, Ei. 
The phenomenon of avoided crossing of resonance states in the complex plane is in complete 
analogy to the well-known phenomenon of avoided crossing of discrete states. In the latter 
case, the crossing point can be found by opening the system and varying the coupling 
strength of the discrete states to the continuum, i.e. by continuing into the complex plane. 
In both cases, the crossing point influences strongly the properties of the states although it 
is hidden 0]. 

The formalism for the description of double QDs with more complicated structure is given 
in We will not repeat it here. We will however use it to obtain some numerical results 
for the transmission through double QDs with a larger number of states. 

III. BRANCH POINTS IN THE COMPLEX PLANE 

Let us define the value 



by which the two eigenvalues ^1^3 differ according to (jH)). F is real only when k = mr; n = 

0. 1. ... When F > 0, Eq. (jH)) gives repulsion of the two levels 1 and 3 in their energies 
Re{zk). When however F < 0, there is a bifurcation of the decay widths Im{zk). 

Most interesting is the case F = since the eigenvalues and eigenvectors of ifefr have 
some special properties under this condition. From (jSJ follows zi = z^ for the eigenvalues, 

1. e. the condition F = defines a point of coalesced eigenvalues. According to Q, the 
components of the (complex) eigenvectors |1) and |3) become infinitely large, and 



Also the normalization condition (lllj) is fulfilled when F = due to the biorthogonality 
of the eigenfunctions, since the difference between two infinitely large numbers may be 
or 1. These relations between the eigenvalues and eigenvectors of Hes that follow from the 
condition F = 0, hold not only for the special case considered here. They hold also for the 




(16) 



l) = ±z|3) 



when F = . 



(17) 
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eigenvectors of an effective Hamiltonian that describes atoms under the influence of a laser 
field 0. Mo.e ge„e.a.., characterize *Ke eigen.ates of a. effec«ve Ha.UtCan that 

npn 

describes an open quantum system j3|, 17|, |20|] . 

The point at which F = 0, is a branch point in the complex plane 

saa. This 

point separates the scenarios with level repulsion on the one hand and width bifurcation 
on the other hand The study on the basis of a schematical model provided the 

following additional results: level repulsion is accompanied by the tendency to reduce the 
differences between the widths of the two states, while width bifurcation is accompanied by 
level clustering. 

According to (jHI), the two eigenvalues zi and of the effective Hamiltonian ((Tj) coalesce 
when Re{F) = and Im{F) = 0. The first condition gives 

vt = i<L^) - elf + (18) 

From the second condition and E = —2cos{k), we find the energy at which the coalescence 
takes place, 

E^^'-mizM, (19) 

In Fig. ^ we present the typical evolution of the real and imaginary parts of the eigen- 
values Zk of the effective Hamiltonian Hcs as a function of the coupling constant v. The 
parameters of the system are e = 2 — L/5, £1 = 1, u = 0.25, L = 3. With these 
parameters, it follows from Eqs. (fTH|) and (fT^ that the eigenvalues zi^s coalesce when 
V = Vc = (1/2 + 9/25)^/^ = 0.9013 and E = = 0.9847. The results shown in Fig. □ 
are obtained for the energy E = E^. Although there are three eigenstates, only zi and 2:3 
coalesce at the point E = Ec, v = Vc- The second eigenstate does not interact with the 
two other ones since it is not directly coupled to the leads. It is coupled to the leads only 
via the two single QDs, and this coupling is symmetrically. This result is in accordance to 
([7j). It can be seen further, that the two states |1) and |3) with energies Re^z^), Re{zi) 
coalesce (when v = Vc) at the energy E = E^. At this branch point in the complex plane 
Ek 7^ Re{zk)\E=E^, El 7^ Re{zi)\E=Ec ■ This means, the two resonance states |1) and |3) 
do cross at E = Ec but not at the energy Ek or Ei. In Fig. |21 (a), the corresponding 
transmission probability versus v and E is shown. 

Let us consider now the behavior of the eigenvalues of the effective Hamiltonian as a 
function of v at the energy E = Ek where Ek = Re{zk{Ek)) is solution of Eq. In the 
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Im(z^) 



FIG. 1: The evolution of Re{zk) (a) and Im{zk) (b), k = 1,3 (solid lines), k = 2 (dashed line), 
of the three eigenvalues of the effective Hamiltonian i/gfr as a function ol v at E = Ec = 0.9847. 
The parameters of the double DQ system are chosen as ei = 1, e(L) = 2 — L/5, u = 1/4, L = 3. 
At V = Vc = 0.9013 the two eigenvalues zi and coalesce. The Re{zi) and Re{z3) approach each 
other when v < Vc, while the corresponding Im{zi) and Im{z^) bifurcate when v > Vc- At the 
branch point in the complex plane Ec^ Ei^., Ei. 



general case, it is not easy to find the solution of the fixed point equation. However for the 
energy at which the eigenvalues coalesce, Eq. ffT^ can be easily solved analytically. 
From ©, (HHI) and (Uni) we obtain 

2{e{L,)-el) 



Ek = e{Lc 



(20) 
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FIG. 2: The transmission probability through the double QD versus v and energy. Each single 
QD has one level at ei = 1. It is e(L) = 2 — L/5 and L = 3. The eigenenergies of the double QD 
are shown by stars. The case (a) corresponds to Fig. ^ The coupling constant between the single 
dots and the wire is u = 1/4. The point of coalesced eigenvalues is Vc = 0.9013, Ec = 0.9847, and 
the solutions of the fixed point equations (|14j) give Ek = Ei ^ Ec as can be seen from Fig. ^ In 
the case (b), the coupling constant u = Uc = 0.1443 between the single dots and the wire is chosen 
in correspondence to Eq. H21|) . Therefore, E^ = Ei coincides with E^ = 7/5 aX, v = Vc- 



With the parameters chosen in Fig. the last equation implies that solutions exist if 
e(L) < 2. We can consider therefore the evolution of the eigenvalues Zk with v at E = Ek = 



and 




(21) 
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e{L) = 7/5 and look for the point where the two eigenvalues coalesce. The critical values 
at the branch point in the complex plane are Uc = 0.1443 and Ec = 7/5. The evolution 
of the eigenvalues Zk with v for Uc, L = 3 is similar to that given in Fig. ^ It is not 
shown here. The corresponding transmission picture Fig. |21 (b) is also similar to Fig. |21 
(a). The main difference is the smaller spreading of the eigenvalues of Hb and the smaller 
transmission probability according to the smaller value u in the case with Ek = Ei = Ec- 
In both cases, the transmission is more spread in energy aX v < Vc than at f > Vc- This 
is in accordance with level repulsion seen in the eigenvalue trajectories at small v and level 
attraction appearing at large v. There is a transmission peak at f ~ 1 near the upper 
border E = 2 in both cases. This peak follows from the energy dependence of the Re{zk): 
the positions of the two resonance states with large width approach E = 2 with f ^ 1 (see 
Fig. ^ where the eigenvalues are shown for an energy E < 2). We can state therefore that 
the characteristic features of the transmission pictures do not depend on whether the two 
states avoid crossing or cross in the complex plane. 

In Fig. ini we present the peculiar symmetrical behavior of the eigenvalues Zk versus v at 
E = for the resonant case with the parameters e(L) = 6i, L = 5. In this case we have, 
according to Eqs. (fTHj) and (fTIJ|) . Ec = and At f < Vc, the widths of the 

two states 1 and 3 are equal, Im{zi) = Im{z-i), while dX v > Vc their positions are equal, 
Re{zi) = Re{z^). The state 2 is not involved in the crossing scenario as in Fig. 

The transmission probability versus energy and v is presented in Fig. IH It has the same 
symmetrical behavior as the eigenvalue pictures. Of special interest is, as Fig. 0] (b) shows, 
that this symmetrical case is at f = 0.53 a perfect filter: the transmission probability is 
equal to one in a large energy range. 

Up to now, we traced the appearance of a branch point in the complex plane by enlarging 
the coupling strength v between system and leads. In such a case, the branch points at which 
two eigenvalues coalesce, appear in a natural manner. It is less evident that the branch points 
in the complex plane can be seen in all parameters of the double QD system that define Eq. 
()18|1 . We can take arbitrary but fixed values of v and u and consider the length L or even 
the energy as a parameter in order to trace the coalescence of zi and z^ at and Ec- 
The corresponding equations for achieving the coalescence are 




(22) 
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FIG. 3: The evolution of Re{zk) (a) and Im(zk) (b), k = 1,3 (solid lines), k = 2 (dashed line), 
of the three eigenvalues of the effective Hamiltonian i^eff as a function ol v at E = Ec = 0- The 
parameters n = 1/4, L = 10, ei = 0, e(L) = 2 — L/5 of the double QD system are chosen in such a 
manner that e(L) = ei = at E = 0. Here, the two eigenvalues coalesce. Vc = S^^^u^^"^ = 0.8409. 

A whole branch cut occurs along L when u = Uc, v = Vc and E = E^. are fixed but is not 
fixed. We consider in the following one branch point corresponding to a fixed value of Ei. 

The case with L as a parameter is shown in Fig. Elfor the same double QD system as 
in Fig. but f = 1. There are two branch points in the complex plane corresponding to 
= V2, Lie = 1-4645 and E2c = -V2, i^2c = 8.5355. When L < Li, and E > V2, the 
two levels 1 and 3 avoid crossing as in the foregoing cases. In the region Lie < L < L2c 
and —a/2 < E < \/2, the levels do not cross at all in the complex plane due to their 
different widths: one of them is trapped by the other one due to the strong interaction via 
the continuum (i.e. via the modes propagating in the leads). For L > L2c and E < —V2, 
the levels again avoid crossing in the complex plane since the widths and with them the 
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FIG. 4: (a) The transmission probability through the double QD versus v and energy for the case 
shown in Fig. |31 (b) The same as (a) but for fixed v = 0.2 (dashed line), v = 0.53 (solid line), and 
V = 0.83 (dot-dashed line). At f = 0.53, the double QD is a perfect filter. 

external coupling of the states via the continuum decrease in approaching E = —2. 

The appearance of two branch points in the complex plane in Fig. El illustrates in a very 
convincing manner the interplay between internal and external interaction in approaching 
a branch point. In any case, a branch point separates regions with avoided level crossing 
{L < Lie, L > L2c) from those without any crossing of the levels {Lie < L < L2c) in the 
complex plane. One should underline, however, that the first branch point influences the 
physical observables such as the transmission probability [Fig. ini(a)], indeed. The second 
branch point occurs as a threshold effect far from the energies Ei and E3 of the two states. 
The two eigenvalues Zi and coalesce at the energy Ec = —V2 <^ E}^ — Vk/2,Ei — F//2, i.e. 
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FIG. 5: The evolution of the real and imaginary parts of the eigenvalues z^, k = 1,3 (solid lines), 
k = 2 (dashed line), as a function of the length L for the same double QD system as in Fig. ^ 
but V = 1. Ec = =b\/2- The critical values of the length L at the two points of coalescence of 
eigenvalues are Lie = 1.4645, Lgc = 8.5355. (a, b) ^ = -^2 - 0.1, (c, d) E = -^/2, (e, f) 
E = -a/2 + 0.1, and (g, h) E = ^/2. 



at the tails of the resonance states. This does not have any influence on the transmission 
probability. 

In Fig. ini(b), the transmission probability is shown at L = L2c- It shows one peak, caused 
by the narrow resonance state, on the background created by the two broad resonance states. 
The narrow resonance is of Fano type by taking into account that the background decreases 
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FIG. 6: (a) The probability T for transmission through the double QD versus E and v for Lc = 
8.5355. (b) The transmission probability as a function of E for fixed v = 0.85. It has one narrow 
peak on the background caused by the two broad resonance states. Parameters: ei = 1, e(L) = 
2 - L/5, u = 1/4 as in Fig. □ 

in approaching the two borders E = ±2. The transmission probabihty for other values of 
L > Lie is similar to that shown in FiglHl 

In Fig. [71 we show the analogue pictures for the E dependence of the eigenvalues z^. 
Due to the fact that the energy is bounded from below {E = —2) and above {E = 2), the 
energy dependence of Im{zk) can not be neglected. It is especially large for states that are 
strongly coupled to the continuum. While the energy dependence of Im{zk) is more or less 
symmetrically around = 0, the Re{zk) show an unsymmetrical behavior as a function 
of energy. It is of special interest, that the branch points in the complex plane appear 
also in the energy dependence of Re{zk) and Im{zk)- An example is the branch point at 
Ec = V2, Lc = 1.4645 that can be seen in Fig. [3 

IV. TRANSMISSION THROUGH A DOUBLE DOT SYSTEM WITH DIFFER- 
ENT COUPLING STRENGTHS TO THE TWO LEADS 

Till now we considered the case that the double QD is coupled to the left and to the 
right reservoir with the same strength v. The couplings may be, however, different from one 
another. Such a case is interesting, also from a theoretical point of view, since the effective 
Hamiltonian becomes unseparable when the two coupling strengths differ from one another. 
This is in contrast to ((Zj) where the double QD is assumed to be coupled symmetrically to 
the reservoirs and, according to (jH)) and ©, the eigenstate |2) does not interfere with the 
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FIG. 7: The evolution of the real (left column) and imaginary (right column) parts of the eigenvalues 
Zk, A; = 1, 3 (solid lines), k = 2 (dashed line), as a function of energy for transmission through the 
same double QD system as in Fig. ^ but u = 1 as in Fig. |SJ The point of coalesced eigenvalues 
is Ec = V2. The critical length is Lc = 1.4645. (a,b) L = - 0.1, (c,d) L = Lc, and (e,f) 
L = Lc + 0.1. 



other two states |1) and |3). 

Following we can write as follows 



\ I ctti / m mn ^ ^^^^ i_2 E + tO - E' 



E 5 



(23) 



where v,w are the coupling strengths between the system and, respectively, the right and 
left reservoirs. Substituting the eigenstates of the closed double QD system © into we 
obtain the following expression for the (symmetrical) effective Hamiltonian 



eff 



)i2 ,,,2 \^ifc 



2T;(r)+Ae) 2y/ri{-ri+A6) 
ui^)e'* / 2 I 2\ ik /n u(v'-'—w'-')e'-' 



\ 



\ 



2^Jr,{r,+Ae) 
2V2ri 



(24) 



2i/r?(»?-A£) 2r,{-n~Ae) j 

The transmission probability for a system with different couplings of the double QD to 
the reservoirs demonstrates new features that appear when v and w differ strongly from 
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one another (Fig. |S)). In the calculations, we have chosen the following parameters for 
the double QD system: e(L) = 2 — L/5, L = 4, u = 0.15, ei = 1. Then from we 
have = 0.8665, = ei = 1, E^ = 1.3345 for the three states of the closed system. 
The positions of the real parts Re{zk),k = 1,2,3, of the three eigenvalues of the effective 
Hamiltonian H^s are given in Fig. |Hl left column, for E = 1.0, 0.92 and 1.26. 

Let us at first tune the energy of the incident particle to be resonant with the eigenenergy 
E = E2 = 1 of the closed system. As it can be seen from Fig. |H1 (a), we can have resonant 
transmission through the system at this energy only for w < 1/2. Correspondingly, the 
transmission probability decreases for large w, Fig. |H1 (b). Next, let us take E = 0.92 that 
approaches Ei for w ^ 1/3 according to Fig. |S1 (c). Resonance transmission through the 
system is possible, at this energy, only when w >l/3 and v = 0.06. Since Re{zi) is almost 
constant as a function of w when w > 1/3, also the transmission remains almost constant 
for w > 1/3. Obviously the transmission is symmetrical relative to v ^ w. As a result 
we obtain the peculiar picture of transmission probability shown in Fig. |H1 (d). A similar 
picture is obtained if the energy is tuned to the third eigenenergy that is E = E3 = 1.26 for 
large w, as shown in Figs. |H1 (e, f). We mention, however, that at larger u the transmission 
picture is less peculiar. Maximum transmission appears when w ^ v and v is about 2 or 3 
times larger than u. 

V. TRANSMISSION THROUGH A DOUBLE DOT SYSTEM WITH MORE 
THAN THREE STATES 

We show now results of some calculations for the transmission through a more realistic 
double QD system with more than one state in each of the single QDs. The number of 
propagating modes in the leads as well as in the wire, connecting the two single QDs, is 
restricted to one as in the foregoing calculations. 

In Fig. ini we show the transmission through such a double QD system with two states 
in each single QD as a function of energy E and length L for u = 0.25 and for four different 
coupling strengths f < 1. The results show the change of the transmission picture as 
a function of L for different v. At small v, the transmission takes place mainly at the 
energies E^ of the discrete states of the double QD. This behavior is called usually resonant 
transmission. At larger v, however, the transmission peaks have nothing in common with the 
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FIG. 8: Left column: The evolution of the real parts of the eigenvalues of ((231) ^ ^ function of w 
for V = 0.1, E = 1.0 (a), v = 0.06, E = 0.92 (c), and v = 0.1, E = 1.26 (e). The parameters 
of the closed double QD system are L = A, u = 0.15, ei = 1, e{L) = 2 — L/5. The circles at the 
x-axes denote the energies E. Right column: The transmission probability through the double QD 
versus coupling v with the left reservoir and w with the right reservoir. The energies E are the 
same as in the corresponding figures of the left column. 

positions Ej^ of the eigenstates of Hb- Here, the energy and L dependence of the transmission 
follows basically that of the wave inside the wire, e = 3/2 — L/7. The transmission picture 
given in Fig. El corresponds to those shown in [l6|. Transmission zeros appear for all v at 
= (el + £|)/2 where ei {k = 1,2; s = 1 , r) are the eigenenergies of, respectively, the 
left and right single QD. It is = = 3/4 in Fig. M 

The eigenvalue pictures corresponding to Fig. El are shown in Fig. [TUl As long as v 
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FIG. 9: The transmission through a double QD versus E and L for v = 0.25 (a), 0.5 (b), 0.75 
(c) and 1.0 (d). The soHd hues represent the five real eigenvalues of the Hamiltonian Hb as a 
function of L. The dashed lines show the eigenenergy of the wire e = 3/2 — L/7. The eigenenergies 
of the two single QDs are equal: £i = 1/2, £2 = and u = 0.25. The transmission zero at 
£"0 = 3/4 is independent of L and v. 

is small, the energies Re{zk) show a dependence on the parameter L that is typical for 
interacting (discrete) states. The Re{zk,i) of the two outermost states avoid crossing at a 
certain L = where the decay widths 2 Im{zk,i) cross. At larger v, however, the eigenvalue 
pictures change since the widths of the two outermost states do no longer cross in the complex 
plane. Though the trajectories projected onto the energy axis cross at a certain value of L, 
the decay widths do not cross at all. This is due to the large difference between Im{zi) and 
Im{z^) as a consequence of resonance trapping (width bifurcation). 

We can see from the eigenvalue trajectories Fig. El that the picture El (d) corresponds 
also to resonant transmission in spite of the fact that its structure is completely different 
from that in El (a). The point is that the eigenvalues of iYefr differ fundamentally from those 
of Hb if the coupling of the states via the continuum is strong. The transmission peak 
appears at the position of a narrow resonance state. Besides this state, there are two broad 
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FIG. 10: The evolution of real (left) and imaginary (right) parts of the five eigenvalues of the 
Hamiltonian H^s as a function of the length L for a double QD system. The coupling of the 
system to the continuum is w = 0.35 (a, b), 0.8 (c, d), and 1.1 (e, f). The parameters of the 
system are u = 0.25, E = 0.25, e = 3/2 — L/7. The energies of the two single QDs are the same: 
El = 1/2, £2 = 1- The transmission of this double QD is shown in Fig. 

and two narrow resonance states lying each very close to one another. The interferences 
between them are obviously destructive. 

Another interesting result seen in Fig. El is that the decay width of the state in the middle 
of the spectrum vanishes at L ^ 3 for all v. At this value of L, the middle state crosses the 
energy -E*^"^ = 3/4 where the transmission is zero. For a discussion of the transmission zeros 
see jl6|. 
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FIG. 11: The transmission through a double QD versus v and E with the length L = 2 (a) and 5 
(b). The parameters are u = 0.25, e(L) = 2 — L/4, £i = 1/2, £2 = 1. The transmission zero at 
£^0 = 3/4 is independent of v and L. 

In Fig. El the transmission through a double QD with altogether five states is shown 
as a function of energy and v for two different lengths of the wire, L = 2 and 5. Each of 
the two single QDs has two levels at ei = 1/2 and 62 = 1, and the mode in the wire is 
e(L) = 2 — L/4. Transmission zeros appear a.t E = 3/4 (for a detailed discussion of the 
transmission zeros see 

The eigenvalue pictures corresponding to Fig. ^2 at E = 0.75 are shown in Fig. E| We 
see a bifurcation of the widths as discussed in Sect. Ill as well as the corresponding branch 
points in the complex plane. At large v, there are two broad resonance states according to 
the two modes propagating in the two leads. The remaining three states are narrow at large 
V. They are trapped by the two broad states. As shown in Fig. El the two outermost states 
coalesce only at L = 3.03. The resonance state in the middle of the spectrum coalesces, 
however, with another state at lower energy for all three lengths L shown in Fig. El 

The eigenvalue pictures calculated at different energies differ from one another in some 
details. The eigenvalue picture [T21 corresponds to Fig. [D calculated at a positive energy 
E. The two broad states are shifted to higher energy when v is large. The shift is in 
the opposite direction when the eigenvalue pictures are calculated at negative energy. The 
calculation at E = gives a symmetrical picture corresponding to Fig. El In this case, the 
positions of all states at large v are almost constant. The resonance trapping mechanism 
occurs symmetrically at -E = 0: the two outermost states coalesce at a somewhat higher 
value of V than the two states lying nearer to the center of the spectrum. The state in the 



23 




0.5 1 1.5 -1 



Re(z,) 



Im(z^) 



FIG. 12: The evolution of real (left) and imaginary (right) parts of the five eigenvalues of the 
Hamiltonian Hcg as a function of the coupling strength v for a double QD. The length of the wire 
is L = 0.7 (a, b), 2 (c, d), and 3.03 (e,f). Parameters: u = 0.25, E = 0.75, e(L) = 2 - L/4, ei = 
1/2, £2 = 1- The transmission of this double QD is shown in Fig. 1111 



middle of the spectrum does not coalesce with any other state. It corresponds to the mode 
moving in the wire and is symmetrically coupled to the states at higher and at lower energy 
when E = 0. This result corresponds completely to those shown in Figs. El 

The figures show clearly that the transmission peaks appear at the positions of the eigen- 
states of Hb only when v is small. At larger v, the transmission is determined by interferences 
between the contributions from the different states. Nevertheless, it is resonant in relation to 
the eigenstates of the effective Hamiltonian Heg. Level repulsion at small v and level attrac- 
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FIG. 13: The transmission through a double QD versus v and E with the parameters L = 1.5 and 
u = 0.2. Each single QD has five levels at = 1/4, 1/3, 1/2, 3/4, 1. The energy in the wire is 
e = 1 — L/8. The four transmission zeros are independent of v and L. 

tion at large v cause features of the transmission pictures for a double QD with altogether 
five states (Figs. ITTlandfT^ that are the same as those of a double QD with altogether only 
three states (Figs. ^tolH). The only difference is the appearance of transmission zeros (Fig. 
ITT|) when the two single QDs are coupled to one another so that the double QD is effectively 
different from a Id-chain as in Figs. ^Jand^l see [l^ . 

In Fig. El the transmission through a QD with five states in each single QD is shown, 
and Fig. El gives the corresponding eigenvalue trajectories of all 11 states. The main 
features discussed for the cases with a smaller number of states remain. This holds true 
also for the transmission zeros the positions of which are determined by the energies of the 
eigenstates of the two single QDs. One of the differences to the cases with altogether three 
or five states is the following. The eigenenergy trajectories aX E = Q are symmetrical around 
the energy = in Fig. (HI with only one state in each single QD, while the symmetry is 
somewhat disturbed in Fig. Elwith more states in each single QD. In the latter case, the two 
outermost states do not approach each other completely. The lower state approaches one of 
the states out of the middle, and the upper state becomes trapped by these two states. As 
a consequence, the region with maximum transmission does not occur in the middle of the 
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FIG. 14: The evolution of the 11 eigenvalues zt of the effective Hamiltonian //off as a function 
oi V aX E = Q. (a) Re{zk), (b) Im{zk). L = 1.5, n = 0.2. Each single QD has five levels at 
Ei = 1/4, 1/3, 1/2, 3/4, 1. The eigenenergy of the wire is e = 1 — L/8. The transmission of this 
double QD is shown in Fig. 1131 

spectrum but at a somewhat lower energy. The reason for this asymmetry is the following: 
the functions Re{zk) of ten states are raising with energy while all the Im{zk) are vanishing 
at the two limits E = ±2 of the energy window (compare Fig. [Zj). Therefore, the widths of 
the states at lower energy are larger than those of the states at higher energy so that they 
trap the higher-lying states. For details of the resonance trapping phenomenon see ItI. 

Common to all the pictures shown in this section is that the single-channel transmission 
through a double QD is of resonant character although its structure depends strongly on the 
strength v by which the dot is coupled to the attached leads. The point is that the evolution 
of the eigenvalues of the effective Hamiltonian ifeff as a function of external parameters 
changes fundamentally at branch points in the complex plane. The transmission through 
the double QD shows a correspondingly sensitive dependence on the external parameters. 
Qualitative changes in the transmission picture are caused by branch points in the complex 
plane which separate the scenario with avoided level crossing from that without any crossing 
in the complex plane. While transmission occurs in the whole energy region with several 
peaks in the case with avoided level crossings, there is a smaller number of peaks of mostly 
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different height in the case without any level crossings in the complex plane. The position 



of these peaks changes as a function of L. Common to both scenarios are only t 



le L 



3 



independent transmission zeros (for a detailled discussion of the transmission zeros see 

The two coupling strengths v and u stand, respectively, for the coupling of the double 
QD as a whole to the leads (environment) and the coupling of the two single QDs to the 
wire (inside the double QD system). The ratio v/u characterizes therefore the ratio between 
external and internal interaction of the states of an open quantum system. When the external 
coupling is much larger than the internal coupling, the external coupling of the levels via 
the modes propagating in the two leads, prevents the formation of a uniform QD. In the 
opposite case of large internal coupling, the relatively weak external coupling is unable to 
break the uniform QD. Most interesting is, of course, the transition region between the two 
different types of bonds in double QDs. 

It is worthwhile to notice the following. The two levels that are the outermost ones of the 
spectrum, cross or avoid crossing in the complex plane aX E = Q. The distance in energy to 
the crossing or avoided crossing, that occurs between two other levels, is smaller than their 
decay widths. That means, effectively all states are involved in the scenario of avoided level 
crossing in the complex plane. 

Additionally, we mention that the dependence of the transmission on the length L of the 
wire is determined by the manner the wave propagates inside the wire. It can be replaced 
by another relation between e and L than that used in our calculations or by the analogue 
relation between e and the width d of the wire. In the last case, L can be kept constant in 
studying the dependence of the transmission from d, see the discussion at the end of Ref. 

Q- 

VI. SUMMARY 

The results considered in the present paper are obtained in the formalism worked out in 
jioj l for the description of a double QD system. The formalism is based on the S matrix 
theory with use of the effective Hamiltonian that describes the spectroscopic properties of 
the open quantum system. The formalism is applied in jl^ to the description of trans- 
mission zeros in the conductance through double QDs. These zeros are determined by the 
spectroscopic properties of the constituents of the double dot system and by the manner 
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the single QDs are coupled. They appear at all ratios v/u of the coupling strengths. Our 
present study is devoted, above all, to the transmission peaks. Their positions and widths 
depend on the ratio v/u and are influenced by branch points in the complex plane. At these 
points, the transition between the two scenarios with avoided level crossing and no crossing 
in the complex plane takes place. In any case, the transmission is resonant. 

As long as v/u is small, the levels repel in energy (as the discrete eigenstates of Hb) and 
the decay widths of the different states are of comparable value. This causes some spreading 
of the transmission probability over a relatively large energy region. At large v/u, however, 
the levels attract in energy and the decay widths bifurcate. This causes transmission peaks 
at the positions of the narrow states that appear on the smooth background created by 
the broad states. The positions of the transmission peaks depend, in this case, strongly on 
the length of the wire or on another parameter that controls the propagation of the mode 
inside the wire. The two different scenarios are separated by a branch point in the complex 
plane. At such a point, two eigenvalues {zk and zi) of the effective Hamiltonian coalesce at 
the energy E = Ec- Sometimes, Ec = E^ = Ei. Mostly however Ek 7^ Re{zk)\E=Ec and 
El 7^ Re{zi)\E=Ec, and the branch point in the complex plane is not a double pole of the S 
matrix. 

We underline that the resonance phenomena appearing in the transmission through dou- 
ble QDs are the same as those observed in, e.g., the scattering on nuclei or atoms 7|. The 
role of the branch points in the complex plane for the transmission through a double dot 
svstem agrees with that discussed in a schematical study and for a double-well system 
In our model double QD, however, the energy dependence of the eigenvalues Zk of 
the effective Hamiltonian H^g is relatively strong. Especially Im{zk) shows a strong energy 
dependence due to the energy window with thresholds at a lower and an upper finite energy. 
The spectrum is therefore bounded from below and from above, and the eigenvalues of the 
effective Hamiltonian cannot satisfyingly be approximated by the poles of the S matrix. 

The results discussed here are true for single-channel transmission through a double QD 
system that consists of two single QDs with similar energy spectra and a narrow wire that 
couples the two single QDs and allows the propagation of only one mode. When the energy 
spectra of the two single QDs are very different from one another and the coupling strength 
u to the wire is small, the transmission picture at large v differs from that discussed above. 
In such a case, the transmission is hindered at large above all due to the energy gap 



28 



between the levels of the two single QDs through which the transmission takes place. 

In the present paper, the behavior of a simple model is considered that reflects many 
characteristic features of realistic double QDs with more complicated structure, see Q|. The 
results obtained may guide the construction of double QDs. The position of transmission 
zeros and transmission peaks can be controlled by varying the coupling strengths v and u 
as well as the propagation of the mode inside the wire. An example is the broad plateau 
with maximal transmission shown in Fig. |l](b). Using the interplay between internal and 
external interaction allows to control the properties of QDs in a systematic manner. 
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